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We study the influence of composition changes on the glass transition of binary hard disc and 
hard sphere mixtures in the framework of mode coupling theory. We derive a general expression 
for the slope of a glass transition line. Applied to the binary mixture in the low concentration 
limits, this new method allows a fast prediction of some properties of the glass transition lines. The 
glass transition diagram we find for binary hard discs strongly resembles the random close packing 
diagram. Compared to 3D from previous studies, the extension of the glass regime due to mixing 
is much more pronounced in 2D where plasticization only sets in at larger size disparities. For 
small size disparities we find a stabilization of the glass phase quadratic in the deviation of the size 
disparity from unity. 

PACS numbers: 64.70.P-, 64.70.Q-, 82.70.Dd 
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I. INTRODUCTION 

Adding a second component to a one-component liquid 
changes its static and dynamical properties. For instance, 
if one adds a low concentration of rather small species, 
depletion forces between the larger particles are induced 
These effective forces are attractive for small sepa- 
rations and tend to stabilize the liquid phase in addition 
to influencing the transport properties @, d, 0). Such 
mixing effects are interesting from a fundamental point 
of view, but also for applications. It is our main goal 
to study the influence of mixing on the glass transition 
of binary systems with hard core interactions in two and 
three dimensions. This will be done in the framework of 
mode coupling theory (MCT) @. 

Mixing effects on the MCT glass transition were stud- 
ied first by Barrat and Latz [6( for binary soft spheres. 
However, the first systematic investigation was per- 
formed by Gotze and Voigtmann [?J for binary hard 
spheres with moderate size ratios 5 = R s /Rb- R s and Rb 
are the radii of the small and big spheres, respectively. 
For size ratios close to unity, a slight extension of the 
glass regime was observed. Larger size disparities induce 
a plasticization effect, leading to a stabilization of the liq- 
uid due to mixing. The results qualitatively agree with 
those from dynamic light scattering experiments [1, 
and molecular dynamics simulations [Tol . [Tl| . In con- 
trast to this, a recent theory of Juarez-Maldonado and 
Mcdina-Noyola based on the self consistent generalized 
Langevin equation (SCGLE) [l2j predicts a plasticiza- 
tion effect also for size ratios close to unity. These au- 
thors argue that the data available from simulations and 
experiments are not sufficiently accurate to rule out one 
of the scenarios. Size ratios far from unity, i.e. 5 <C 1, 
may be problematic. First, the quality of e.g. Percus- 
Yevick (PY) theory used to calculate the static input for 
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MCT may become less reliable. Second, phase separa- 
tion (see discussion in Ref. [l3[) and third, a sequential 
arrest of the big and small particles (by a type-A transi- 
tion) could occur. The diverging lengthscale associated 
with a type-A transition affects the quality of the MCT 
approximations. 

The results of Gotze and Voigtmann Q exhibit four 
mixing effects, two of which were mentioned above. 
The two remaining mixing effects are an increase of 
the plateau values of the normalized correlation func- 
tions for intermediate times for almost all wave numbers 
upon increasing the concentration of the smaller parti- 
cles, and a slowing down of the initial part of the relax- 
ation of the big-big correlators towards these plateaus. 
Our motivation is twofold. First, we want to explore 
whether these effects also exist in a corresponding two- 
dimensional liquid of binary hard discs. A recent ex- 
periment [14j has given evidence for glassy behavior in a 
similar two-dimensional liquid including dipolar interac- 
tions. Second, we will investigate in more detail the in- 
fluence of mixing close to the monodisperse system, i.e. 
fixing the packing fraction ip at ipf, (the critical packing 
fraction of the monodisperse system), how does a very 
small perturbation of the monodisperse system influence 
the glass transition? The arbitrary small perturbation 
can be achieved in three ways, either by adding a very 
small concentration of smaller or bigger species for given 
arbitrary S < 1 , or by a slight decrease of the diameter of 
an arbitrary concentration x s of the smaller particles, ac- 
companied by a slight increase of the remaining particles, 
i.e. 1 - S < 1. 

The mixing effects in the low concentration limits fol- 
low directly from the slopes d<p c (x s ,5)/dx s at x s = 
and x s = 1 of the glass transition lines (GTLs) (p c (x s , S) 
at fixed 5. If dip c (x s , 5)/dx s \ Xs =o is positive (negative), 
the liquid (glass) is stabilized. The same is true if 
dip c {x s ,5)/dx s \ Xa= i = -d(p c {l - x b ,5)/dx b \ Xb = is neg- 
ative (positive). Xb = 1 — x s is the big particle's con- 
centration. Since the determination of these slopes from 
the numerical result for tp c (x s ,8) with discretized values 
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of x s is not precise, particularly for 5 closer to unity (cf. 
the critical lines for 5 = 0.7 and 0.8 in Fig. 1 of Ref. 0]), 
we will derive an analytical expression for d(f c (x s , 5) / dx s 
for arbitrary x s and 6. Applied to x s — and x s = 1, 
only the glass transition singularity of the monodisperse 
system is needed. The remaining quantities entering the 
slope at x s = and x s = 1 can be determined from a 
perturbational approach discussed below. The applica- 
tion for the slope formula will be done for both, hard 
discs and hard spheres. This allows to explore the di- 
mensional dependence (at least for d — 2 and d = 3) of 
the mixing effects in the weak mixing limit. 



II. MODE COUPLING THEORY 

We will restrict ourselves to the essential equations to 
keep our presentation self-contained. For details, the 
reader may consult Ref. 5]. Correlation functions are 
matrix valued vectors denoted by bold symbols A, B etc. 
Their components A k , B k being M x M matrices (A? ), 
(B^) (in case of an M-component fluid) are labelled by 
subscript Latin indices (the wave numbers) which can be 
taken from a discrete or a continuous set. The elements 
Afc, B^f of these matrices are indicated by superscript 
Greek indices, in some cases these elements shall also be 
denoted by (A)^ , (B)^f . Matrix products are defined 
componentwise, i.e. C = AB reads C k = A^B^ for all 
k. We call A positive (semi-) definite, (A>zO), A >- if 
this is true for all A k . denotes the (generalized) zero 
matrix. If k is restricted to a finite number of values, 
then the standard scalar product of A and B shall be 
defined as (A\B) = J2 k J2 a , Af Bf . 



A. General equations 



5af3 / (k 2 D^x a ) where hydrodynamic interactions are ne- 
glected. £)° denotes the short-time diffusion coefficient 
of a single particle of the species a inserted into the fluid. 
With this, the short-time asymptote of $(t) is given by 



*(* -> o) = s - r~H + o{t 2 ), 



(2) 



here we restrict ourselves to t > 0. MCT approximates 
m(t) by a symmetric bilinear functional T of $(t), 



m(t)=^[#(*), #(*)]. 



(3) 



It is straightforward to generalize the explicit expression 
for T of a simple fluid in d dimensions (d > 2) presented 
in Ref. [151 ] to multicomponent systems. The result is 
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where 



_ (k 2 +P 2 - q 2 )cfS al + (fc 2 -p 2 + q 2 )c^8 a0 



kpq 



[4k 2 p 2 - (k 2 + P 2 -q2)2](3-d)/4 



(6) 



c^f denote the direct correlation functions, c is related 
to S via the Ornstein-Zernike (OZ) equation 

a/3 



(S = S afj /x a -nc%- 



(7) 



n is the total number of particles per volume and Qd = 
2-K d / 2 /Y(d/2) the well known result for the surface of a 
unit sphere in d dimensions. T(x) is the gamma function. 



We consider an isotropic and homogeneous classical 
fluid consisting of M macroscopic components in d di- 
mensions. <&(i) denotes the matrix of time dependent 
partial autocorrelation functions of density fluctuations, 
4>£^(t) {a, (3 — 1,...,M), at wave number k. We re- 
quire the normalization 4>(0) = S, where S denotes 
the static structure factor matrix whose elements obey 
linifc^oo = x a 5 a p. Here S a p denotes the Kronecker 
delta and x a the particle number concentration of com- 
ponent a. 

Considering overdamped colloidal dynamics, the 
Zwanzig-Mori projection operator formalism yields the 
equation of motion 



r*(t) + S- 1 <f>{t) + / dt'm(t - 







(1) 



with the memory kernel m(t) describing fluctuating 
stresses and playing the role of generalized friction, r 
is a positive definite matrix of microscopic relaxation 
times. Its components shall be approximated by = 



B. Definition of the model 

The M-component MCT in d dimensions shall be ap- 
plied to binary hard "sphere" mixtures (HSM) in d di- 
mensions consisting of big (a = b) and small (a — s) par- 
ticles. Let R a denote the radius of the species a. Three 
independent control parameters are necessary to char- 
acterize the thermodynamic state of a HSM. We choose 
them to be the total packing fraction tp = (f s + ifb with 
(fi a = nXafod/d^Ra, the size ratio 6 — R s /Rb < 1, and 
the particle number concentration x s of the smaller par- 
ticles. 

For the following, we discretize the MCT equations, 
i.e. k is discretized to a finite, equally spaced grid of 
K points, k = (pa + k)Ak with k = 0, 1, . . . , K — 1 and 
< da < 1. The integrals in Eq. dU) are then replaced 
by Riemann sums 



/■k+p 
dq . 
,k-p\ 



K—l roin{K— l,fe+^} 

(a^) 2 E E ■■■ ( 8 ) 

33=0 q=\k-p\ 
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and Eq. ([T]) represents a finite number of coupled nonlin- 
ear "integro" -differential equations. We further restrict 
our numerical studies to the cases d = 2 and d = 3. For 
the offset, following previous works, we choose 62 = 0.303 
for d = 2 [IH and 63 = 0.5 for d = 3 Q. The choice 
K = 250 and Afc = 0.3 turns out to be sufficiently accu- 
rate to avoid larger discretization effects. 

For calculations with finite concentrations of both par- 
ticle species, the unit length shall be given by the di- 
ameter 2Rb of the bigger particles, and the short-time 
diffusion coefficients D° shall be assumed to obey the 
Stokes-Einstein law. Further, the unit of time is chosen 
such that D° = 0.01/ (2R a ). For the numerical solution 
of Eq. (JTJ) we use the algorithm first published in [lj|. 
Our time grids consist of 256 points, as initial step size 
we choose 10 -8 time units. 

For the discussion of the weak mixing limits (see be- 
low) it is convenient to choose the diameters 2R a of the 
majority particle species as unit length. 



C. Static structure 

Approximate closures of the OZ equation provide the 
most powerful methods currently available for a fast cal- 
culation of the pair correlation functions from first prin- 
ciples (l7j |. The OZ equation for an arbitrary mixture is 
given by 

h = c + ncxh (9) 

where x^f = x a 5 a p and the h^P are the total correlation 
functions. For our binary HSM model we use the PY 
approximation given by 

h a P{r) = -1, r<(R a + R p ), 

c?*f>(r) = 0, r > (R a + R ). [W > 

In odd dimensions the coupled Eqs. © and (|10[) can be 
solved analytically [lg|. In even dimensions numerical 
methods must be employed. Among the several exis ting 
algorithms [l9j we use the classical Lado algorithm |20l ] 
for simplicity. In our numerical solution of the 2D system 
we use a real space cutoff r max — 50 with 4000 grid 
points. 



D. Glass transition lines 

The nonergodicity parameters (NEPs) F = (F?^) are 
given by F — lim t ^oo <fr(t). For the discretized model 
described above, the following statements can be proved 
[2l| . Equation |T]) has a unique solution. It is defined for 
all t > and is completely monotone, (—d/dt) n ^(t) >z 
0. F y is (with respect to >z) the maximum real, 
symmetric fixed point of the nonlinear map 

X[X] = s -(s- 1 +F[X,X})-\ (11) 



Iterating Eq. (fTTj) starting with X = S leads to a mono- 
tonically decaying sequence converging towards F. Lin- 
earization of X yields the positive definite linear map 
(stability matrix) 

C[Y] = 2(S-F)F[F,Y](S-F) (12) 

with C[Y] y for all Y >z 0. From a physical point 
of view, it is reasonable to assume that C is irreducible 
if F >- 21]. C has then a nondegenerate maximum 
eigenvalue < r < 1 with a corresponding (right) eigen- 
vector H >- 0. For any other eigenvalue f of C, \f \ < \r\ 
holds, and if |f| = |r|, then the corresponding eigenvector 
can not be positive definite. Hence, possible MCT sin- 
gularities are identified by r = 1 and belong to the class 
A t , 1 = 2,3,..., introduced by ArnoPd [13]. The adjoint 
map C of C satisfies (C[A]\B) = [A\C[B]) for all A, 
B . Its eigenvector H >- is the left eigenvector of C cor- 
responding to the eigenvalue r. These two eigenvectors 
are determined uniquely by requiring the normalization 

(H\H) = (H\H{S - Fj^H) = 1. (13) 

For binary HSM models, higher order singularities may 
occur for large size disparities where the packing contri- 
butions i a = (f a ftp of both components are of the same 
order [23]. In the present paper, we restrict our discus- 
sion to the generic (type-B) MCT bifurcations belonging 
to the class A2 where F jumps from to F c >~ 0. Quanti- 
ties taken at critical points shall be indicated by a super- 
script c. The glass transition takes place at the critical 
surface tp c (x s ,S) within the three-dimensional physical 
parameter space (ip, x s , 5) . ip c (x s ,S) fulfills 

<p c (x s ,6) = <p c (l-x s ,l/5). (14) 

Equation (fT4| demonstrates that (p c (x s ,S) for fixed S is 
not symmetric with respect to the equimolar concentra- 
tion x s — 1/2. However, for x s — 1/2 and small dispar- 
ity, i.e. e = (1 - 5) <C 1, it follows from Eq. (Jl4j) that 
ip c (l/2, 5) ~ e 2 in leading order in e. Accordingly, for the 
equimolar situation and small disparity the influence of 
disparity is quadratic only. ip c (x s ,5) can be determined 
numerically by a simple bisection algorithm monitoring 
the NEPs. 



E. Slope of a critical line 

For a general model system with L external, i.e. phys- 
ical control parameters £ = . . . , the generic glass 
transition singularities form a (L — l)-dimensional hy- 
persurface H. Locally, this surface can be represented, 
e.g. as £f(£i, . . . . . . ,£ L ) for any I. For fixed 

i 3, i I, £f(£ii •••)&)•• ■ ,6-1,6+1, •••,££,) de- 
scribes a GTL which is a function of An expression 

for its slope (d£f/<9*y)(£i, 6-1,6+1, ■ ■ ■ , £i) is 

obtained by use of the separation parameter a. Let 
£ c € H be a critical point and A£ = £ — £ c . Then the 
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separation parameter is a linear function <x(A£) in A£*[j|. 
<t(A£) = defines the tangent plane of the hypersurface 
TL at the critical point £ c . Then it is easy to prove that 



• si 



da/d{AQ 



,8) 



(15) 



A£=0 



The separation parameter <r(A£) follows from 



CT(<f) = (ff c |{S c -F c }S c_1 {S^[FSF c ](S-.F c ) 

-S C ^ C [F C ,F C ](S C -F C )}) (16) 

by expanding around £ c up to linear order in A£ 5, 24]. 
The result (fT5|) demonstrates that the separation param- 
eter besides being a measure for the distance from the 
critical point £ c also contains local information of a GTL. 
Applied to a binary liquid, Eq. (p~5j) yields 



dip c 
dx s 



da/d(Ax s ) 



da/d{A<p) 



(17) 



(Aip,Ax B ,AS)=0 



where AS and all critical input parameters have to 
be considered as fixed constants when calculating the 
partial derivatives. A similar expression follows for 
(dip c /dS)(x s = x%5 = 5 c ). 

Let us further remark that the concept of introducing 
a separation parameter is not restricted to MCT mod- 
els. Thus, Eq. (fT5"]) holds for any system which has at 
least two control parameters and exhibits the generic A2 
bifurcation scenario 1 221 . 



F. Weak mixing limit 

One of the central aspects of our paper is to demon- 
strate the predictive power of Eq. (jTTJ) for the limits 
x s = and Xb = 0. By performing these limits ana- 
lytically, we obtain formulae whose numerical evaluation 
is much less time consuming then the numerical proce- 
dure mentioned above, i.e. to determine the slope from 
(p c (x s , 5). Note, the knowledge of the initial slopes of the 
GTLs for both limits is already sufficient to estimate their 
qualitative behavior under certain assumptions. The es- 
sential steps for the calculation of the slope dip c /dx s are 
explained in Appendix [A"l 



III. RESULTS AND DISCUSSION 

A. Glass transition lines 

Figure [T] shows normalized slopes of the GTLs at 
x s = as functions of 1/5 for the binary 2D and 3D 
HSM models. Because 5 = 1 represents a one component 
system, the slopes have to be zero at this point. While 
the numerical results for the 3D model clearly support 




FIG. 1: (Color online) Normalized slopes of the GTLs at 
x s = for the binary HSM models in 2D and 3D. It is ipo — 
0.6914 for d = 2 and 0.5159 for d = 3. For 1/5 > 3.5 in the 
3D model, the tagged particle NEPs indicate a derealization 
transition of the smaller spheres. This regime is indicated by 
open symbols (see text). 



0.09 - 



B 0.06- 



0.03 



-0.03 - 




FIG. 2: (Color online) Slopes of the GTLs at x b = 0. 



this statement, the numerical data for the 2D model at 
5 = 1 slightly deviate from zero (see also Fig. [2]). This, 
however, is an artifact due to the numerically calculated 
static structure factors in 2D. For the 3D model, we have 
used the analytical solution of Eqs. and (jTUJ) to calcu- 
late the static input for MCT which has led to a better 
self-consistency at 5 = 1 then for the 2D model. For 5 
close to unity, the slopes become negative which means 
that the presence of a small concentration of the smaller 
particles stabilizes the glass. After exhibiting a minimum 
at 5—, the slopes become zero again at 5q and remain 
positive for 5 < 5q. Here the presence of the smaller par- 
ticles stabilizes the liquid which is nothing but the well 
known plasticization effect. Upon further decreasing of 
5, the slopes exhibit a maximum at 5+ and indicate a 
monotonic decay for asymptotically small 5. For the 2D 
model, this decay is more stretched then for the 3D case. 
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FIG. 3: (Color online) (a) Relative variation of the glass tran- 
sition lines for the binary HSM model in 2D. The open squares 
calculated with K — 400 grid points (instead of K = 250) give 
an estimate for the error due to the high wave number cut- 
off, (b) Relative variation of the glass transition lines for the 
binary HSM model in 3D using the numerical data of Gotze 
and Voigtmann [7| calculated with Ak — 0.4 and K — 200. 



For the 3D model, we observe a continuous transi- 
tion of the tagged particle NEPs (F c '^)l s (see Appendix 
[Mil to zero by approaching i5 w 1/4 from above. This 
indicates a derealization transition of the smaller spheres 
in the glass formed by the bigger ones [1^, [2(| [53] ■ Such 
a transition is strongly influenced by a l/fc 2 -divergence 
of the memory kernel for the tagged-particle correlators 
at k = [28[ . This singularity reflects the fact that, in- 
side a fluid, the momentum of a single tagged-particle is 
not conserved. Although the evaluation of Eq. (|17p at 
x s = requires {F c ^) s ^ as input, the qualitative in- 
dependence of ip c should not be influenced by this prob- 
lem. Nevertheless, we show the corresponding data for 
1/5 > 3.5 in Fig. Q] with open symbols. However, these 
data show the same qualitative behavior as the corre- 
sponding ones for the 2D model. For our choice of the 
lower cutoff for k, the MCT model does not yield a delo- 
calization transition in 2D, even if we use the PY result 
for R s — as static input. This, however, is an artifact 
due to the singularity of the tagged-particle memory ker- 
nel at k = 0. Again, the qualitative ^-dependence of ip c 
should not be influenced. 

Figure fj] shows normalized slopes of the GTLs at 
Xb = 0. For 8 close to unity, the presence of a small con- 
centration of the bigger particles leads to a stabilization 
of the glass. The slope vanishes at <5q > So. A strongly 
increasing plasticization effect occurs for smaller 8. 

Apart from the problems discussed above, the results 
shown in Figs. [1] and [2] allow to predict the shape of the 
GTLs. Both x s = and Xb = define one component 
models with the same critical packing fraction ipQ . Hence 
the GTLs show a single minimum for S' < 8 < 1, ex- 
hibit a minimum followed by a maximum (S-shape) for 



for smaller 8 < 8q. Here we assumed that two or more 
minima (maxima) do not occur. Figure [3] (a) shows the 
relative variation (ip c — ip^) / 1 ip^ of the GTLs for the binary 
2D HSM model. Results of Gotze and Voigtmann for the 
3D model [7| are shown in Fig.[3](b). The ^-dependence 
of these GTLs agrees with the 5-dependence predicted 
from the slopes at x s — and Xb = 0. Particularly, the 
S-shapes of the GTLs for So < 8 < S' Q are reproduced. 

All our results predict the following trend: Compared 
to the 3D model, the stabilization of the glass is much 
more pronounced in the 2D model where the less pro- 
nounced plasticization effect only sets in at larger size 
ratios. The maximum relative decrease of ip c occurring 
at S = 0.7 in 2D (see Fig. [3] (a)) is about five times larger 
then the maximum downshift of (p c in 3D which occurs 
at 5 = 0.8 (see Fig.[3](b) and Fig. 2 in Ref. Q). Qualita- 
tively, the binary hard disc liquid exhibits the same two 
mixing effects discussed in the introduction for the hard 
sphere liquid. 

A finer resolution of the slope s c (x s ,5) = 
(d(f c /dx s )(x s ,5) in Fig. Q] for delta close to unity 
(see inset) shows that (ds c / d8)(x s = 0, 8 = 1) = 0. The 
numerical data for the 3D model show this behavior more 
clearly then the corresponding ones for the 2D model, 
for technical reasons mentioned above. The resolution of 
Fig. [U already exhibits that (ds c /dS)(x s = 1,6 = 1) = 0. 
Therefore, ip c (x s , S) at x s — and x s — 1 is quadratic in 
e = (1 — 5) for 8 close to unity. Since Eq. (fLij) has led to 
the same £-dependence at x s = 1/2, we conjecture that 



<p c {x s ,8)~{l-Sf 



(18) 



for all x s and small size disparity. A numerical check for, 
e.g. x s — 1/4, has confirmed the validity of Eq. (fTFj) for 
d — 2 and d = 3. Equations (14]) and (TH]) imply 



ip c (x s ,S) = cp c (l - x s ,S). 



(19) 



intermediate 6q < 6 < 6' , and show a single maximum 



Consequently, the GTLs become symmetric in x s with 
respect to x s = 1/2 in the limit of small size dispar- 
ity. Then the maximum enhancement of glass formation 
occurs at equimolar concentration x s = 1/2, excluding 
again the occurrence of more then one minimum. 

Okubo and Odagaki [2^] have numerically calculated 
random close packing values [4l[ ip rcp of binary hard 
discs by use of a so-called infinitesimal gravity proto- 
col. Figure 2] presents their results for (ip rcp — <p r cp ) / r o cp . 
0rcp ^ 0.8139 is close but not identical to the averaged 
value (fiQ Cp = 0.82 for monodisperse hard discs. Despite 
the large numerical uncertainty at x s = 0, x s — 1 and 
5=1 (this might result from the fact that for monodis- 
perse hard discs the applied procedure tends to build up 
locally odered structures), the data show a striking sim- 
ilarity to {ip c — tpg)/ipg (Fig. [3] (a)). The change from 
the single minimum shape to an S-shape and a maxi- 
mum shape by decreasing 5 is clearly reproduced by the 
random close packing result. 
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mensions, we have 



FIG. 4: (Color online) Relative variation of the random close 
packing fraction using the numerical data of Okubo and Oda- 
gaki [291 ] . </?g Cp is determined such that the relative variation 
vanishes below but close to x s = 1. 
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FIG. 5: (Color online) (a) Glass transition lines for the bi- 
nary HSM model in 2D, plotted as functions of the packing 
contribution of the smaller particles x a = <Ps/<fi- (b)-(c) ot- 
relaxation times defined by $ b k b ((T rel ) b k b ) = 0.1(F c ) b k b for the 
correlators of the big particles at k = 5.1909 for fixed 5 and 
p close below the corresponding GTLs for the binary HSM 
model in 2D. In both (a) and (c), the open triangles calcu- 
lated with K — 400 grid points (instead of K = 250) give an 
estimate for the error due to the high wave number cutoff. 



B. Mixing scenarios 



In this section we will demonstrate that the mixing 
scenarios presented in Ref. [7( for binary hard spheres 
are also observable for binary hard discs. For this pur- 
pose, we follow Gotze and Voigtmann [3] and choose <p, 
5, and the packing contribution of the smaller particles 
x s = fs/f a $ independent control parameters. In d di- 



1 + x s (l/5 d - 1)' 



(20) 



As a direct analogon to Fig. 1 in Ref. 0, Fig. [5] (a) shows 
GTLs for the binary HSM model in 2D, plotted as func- 
tions of i s for three representative values for 5. The GTL 
for (5 = 5/7 shows a single, clearly pronounced minimum, 
the line for 5 = 1/2 is S-shaped, and the GTL for 5=1/3 
exhibits a single maximum. 

For both the hard sphere and the hard disc system the 
relative variation of tp c with concentration is of the order 
of one percent or less (see Figs. [3] and [5] (a)). This can 
neither be observed by experiments nor by simulations. 
As already stressed in Ref. p}, the variation of ip c with, 
e.g. x s , may be reflected by a strong variation of the 
ct-relaxation time {T re[ )^f , i.e. a variation of the charac- 
teristic time scale for the final decay of & k ^ (t) to zero in 
the liquid phase. If ip is fixed below but sufficiently close 
to ip c (x si 5), i.e. if <p is fixed such that there exists an in- 
terval in the (x s , <5)-plane such that < <p c (x s ,5) — p <C 1 
is satisfied for all (x s , 5) within that interval, then the ct- 
relaxation time (T rel )^f ~ (tp c (x s , S) — ip)~i( Xs ^ is ex- 
tremely sensitive to the variation of (x s , 5) within that 
interval. Figure [5] (b) shows a-relaxation times defined 
by $fc 6 ((r re 0fe b ) = 0.1(F c )£ b for the unnormalized corre- 
lators of the big particles at k = 5.1909 for fixed S = 5/7, 
1/2 and fixed ip = 0.686 below but close to the corre- 
sponding GTLs for the binary HSM model in 2D for 
different packing contributions x s . The qualitative in- 
dependencies of the corresponding GTLs in Fig. [5] (a) 
are clearly reflected by the i s -dependencies of the a- 
relaxation times. (T rel ) b k b shows a single maximum for 
6 = 5/7, and is S-shaped for 5 = 1/2. For 5 = 5/7, 
{T rel ) h k h varies by more then three decades. Figure [5] (c) 
shows (T rel ) b k b at k = 5.1909 for fixed 5 = 1/3 and fixed 
<p = 0.691 below but close to the corresponding GTL 
for the binary HSM model in 2D. The qualitative in- 
dependence of the corresponding GTL in Fig. [5] (a) is 
reflected by a single minimum in (T rel ) bb . Note that for 
this 5 we had to choose a slightly larger value for ip than 
for the two other examples shown in Fig. [5] (b) in order 
to clearly observe this effect. In contrast to this, Fig. 11 
in Ref. [7[ exhibits all three scenarios for one common <p. 
In our 2D model, however, the minimum of <p c occurring 
for 5 = 5/7 is more strongly pronounced then the corre- 
sponding one for 5 = 0.8 in 3D shown in Fig. 1 in Ref. 0]. 
This fact makes the choice of a common <p for all three 
considered values of 5 for the 2D model difficult. 

Let us also discuss some representative correlators in 
more detail. As a direct analogon to the upper panel 
in Fig. 8 in Ref. , Fig. [6] shows normalized correlators 
® h k h {t)/Sf of the big particles for the binary HSM model 
in 2D at fixed p = 0.686, 5 = 5/7 and k = 5.1909 for 
different packing contributions x s of the smaller discs. 
Let (f rel ) b k b be the characteristic time scale specified by 
90% of the decay from the normalized plateau value 
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FIG. 6: (Color online) Normalized correlators of the big par- 
ticles for the binary HSM model in 2D at tp = 0.686, 5 = 5/7 
and k = 5.1909 for different packing contributions x s = p s /<fi- 
Filled diamonds mark the crossings of the normalized criti- 
cal plateau values (-F e )| / '{S c y^ . Open diamonds mark the 
crossings of the values 0.1(F c ) b k b /(S c )f . 
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FIG. 7: (Color online) The same as Fig. ©but for 5 = 1/3 and 
= 0.691. The thin lines show the short-time asymptotes 
given by Eq. ([2} for x s = 0.01 and x s = 0.3 (from left to 
right). 



{F c )f/(S c )f to zero. For the chosen value of S the 
corresponding GTL shows a single minimum shape, see 
Fig. 03(a). Hence, starting from the almost monodisperse 
system at x s = 0.01 and increasing the packing contribu- 
tion of the smaller discs to x s = 0.3 leads to a decrease 
of the distance ip c (x s , S) — ip from the GTL. This fact is 
reflected by an increase of {T rel ) hb by more then three 
decades (see the open diamonds in Fig. 03) ■ 

An analogous scenario to the upper panel of Fig. 9 in 
Ref. @ is presented in Fig. [7] It shows normalized corre- 
lators $> b k b (t)/Sl b of the big particles for the binary HSM 
model in 2D at fixed <p = 0.691, 5 = 1/3 and k = 5.1909 
for different packing contributions of the smaller discs. 
For the 8 chosen here the corresponding GTL shows a 



FIG. 8: (Color online) Normalized critical NEPs for the big 
particles at 5 = 1/2 for the binary HSM model in 2D. 

single maximum shape (see Fig. 03 ( a ))- Hence, starting 
at x s = 0.01 and increasing the packing contribution of 
the smaller discs to x s = 0.3 leads to an increase of the 
distance ip c (x s , 8) — <p from the GTL. As a result, {f rel ) h k 
decreases by about two decades (see the open diamonds 
in Fig. 13). 

Two additional mixing effects (briefly mentioned in the 
introduction) were reported in Ref. [7j for the 3D model. 
The first of these effects is the increase of the normal- 
ized critical NEPs (Debye- Waller factors) {F c )% a / {S c )% a 
upon increasing x s for almost all k (related to an increase 
of the plateau values of the correlation functions for in- 
termediate times). The origin of this effect is explained 
in great detail in Ref. Q. The 2D model shows simi- 
lar behavior. Here we restrict ourselves to a representa- 
tive example. Figure 03 shows normalized critical NEPs 
{F c )f/{S c )f for the big particles at 8 = 1/2 for the bi- 
nary HSM model in 2D. The data for x s = represent 
the Debye- Waller factors of a monodisperse system con- 
sisting of discs with diameter one, while the result for 
x s = 1 corresponds to the critical tagged particle NEPs 
(Lamb-Mofibauer factors) of a single disc of diameter one 
inserted into a monodisperse system consisting of discs 
with diameter 8 < 1. These Lamb-Mofibauer factors for 
all k are larger then the corresponding Debye- Waller fac- 
tors of the system of monodisperse discs with diameter 
one. Provided that {F c )f/{S c )f varies smoothly for all 
< x s < 1 (i.e. there are no multiple glassy states for 
the considered value of 8), one obtains an increase of the 
Debye- Waller factors upon increasing x s as an overall 
trend (see also the filled diamonds in Figs. [S] and [7]) . 

The second remaining mixing effect is the slowing down 
of the initial part of the relaxation towards the plateau 
values for the correlators of the big particles in the sense 
that 3>| 6 (t)/<Sfc 6 versus log 10 (t) becomes flatter upon in- 
creasing x s . This effect is clearly visible is Figs. 03 and [7] 
The authors of Ref. [3] conclude that the change of the 
short-time dynamics upon increasing x s is not sufficient 
to explain the observed effect. Figure supports this 
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statement. The shown short-time asymptotes resulting 
from Eq. @ for x s — 0.01 and x s = 0.3 fall already at 
log 10 (f) w — 1 significantly below the corresponding cor- 
relators. Thus, the enormous flattening of the curves in 
the region < log 10 (t) < 2 can not be simply explained 
by the slowing down of the diffusion at short times. 

Let us conclude at this point that we have found the 
same four mixing effects for binary hard discs as have 
been reported for binary hard spheres in Ref. 0. The 
subtle scenario in Fig. [7] is the result of an interplay of 
three of these mixing effects. The increase of x s leads first 
to both an increase of the plateau values of the correlators 
at intermediate times and a slowing down of the initial 
part of the decay toward these plateaus. However, the 
increase of x s also leads to a decrease of the a-relaxation 
times, i.e. an enhancement of the final decay to zero, and 
thus to a crossing of the correlators. 



IV. SUMMARY AND CONCLUSIONS 

In the present paper we have studied the influence of 
composition changes on the glass transition for binary 
hard disc and hard sphere mixtures in the framework of 
MCT. 

By deriving Eq. (fl~5|) , we have shown that the well- 
known separation parameter not only describes the scal- 
ing of the NEPs in the glass 0, [3(3] , but also describes the 
local variation of the GTLs to linear order. For low con- 
centration limits of one particle species we have evaluated 
the slopes of the GTLs, Eq. ifTTj) . by using a perturba- 
tion ansatz. With this we have introduced a new method 
which allows a fast prediction of some qualitative proper- 
ties of the GTLs. Note that this method can be applied 
to any MCT model with more then one control param- 
eter. For instance, a similar analysis should be possible 
for hard spheres with attractive potentials in the limit 
of vanishing attraction strength [3l[ , or equivalently, for 
temperature going to infinity. More generally, Eq. (|15p 
holds for any system which has at least two control pa- 
rameters and exhibits the generic A2 bifurcation scenario 

The direct comparison of the models in 2D and 3D 
show similar qualitative behavior. Particularly, the same 
four mixing effects have been found as for hard spheres 
Q. However, we have also found some differences. The 
main difference is the fact that the extension of the glass 
regime due to mixing for size ratios close to unity is more 
strongly pronounced in 2D then in 3D. 

For small size disparity we have presented analyti- 
cal and numerical evidence that the stabilization of the 
glassy state is quadratic in (1 — S) and that the GTLs are 
almost symmetric with respect to their equimolar con- 
centration x s = 1/2. At this concentration the stabiliza- 
tion is maximal. These properties have not been noticed 
before. 

Finally, we have shown that the qualitative (x s , in- 
dependence of ip c for some representative values of S is 



identical to that of the random close packing ip rcp . This 
is particularly true for the S-shape dependence for in- 
termediate values for 5. The maximum shape variation 
of ip c which implies stabilization of the liquid state and 
which has been related to entropic forces 0, HH, [33| 
exists also for <p rcp for smaller S. Since the random close 
packing procedure of Ref. [2^ | is a nonequilibrium pro- 
cess which maximizes the density locally, it is not obvious 
that the stabilization effect is of entropic origin, at least 
for S not too small. 

At this point we should also remember that ip rcp is 
not uniquely defined. For instance, a subsequent shak- 
ing of the configurations produced by the infinitesimal 
gravity protocol used in Ref. [29J would typically lead to 
random structures at even higher densities. Hence, one 
may ask whether the qualitative trends shown in Fig. [4] 
are reproducible by using different procedures for calcu- 
lating ip rcp . A different approach is the investigation of 
jamming transitions of hard discs or hard spheres. Simu- 
lations on frictionless systems of repulsive spherical par- 
ticles have given evidence for a sharp discontinuity of 
the mean contact number Z at a critical volume fraction 
tpjam H,lll. These results are supported by ex- 

periments on binary photoelastic discs with 5 ks 0.86 and 
x s = 0.8 [38[ . Recently, tp^ am has been determined by 
Stark, Luding, and Sperl as function of x s for different 
values for <5, both by experiments on photoelastic discs 
and by corresponding computer simulations [3^ |. Their 
results clearly support all the qualitative features pre- 
sented in Fig. 21 whereby supporting the results shown 
in Fig. [3 

Let us conclude with some open questions which are 
worth to be investigated in the future. For the 3D model, 
higher order singularities (connected to the existence of 
multiple glassy sates) occur below S w 0.4 [23|]. The ques- 
tion, whether such transitions also exist in 2D, requires 
a more detailed numerical study. The consistency of our 
MCT results with the corresponding random close pack- 
ing data supports the quality of MCT in 2D. However, 
also a quantitative comparison of the dynamical MCT re- 
sults with molecular dynamics simulations is necessary. 
A further step towards reality will be the study of MCT 
for binary discs including dipolar interactions for which 
detailed experimental studies exist [Til ]. 
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APPENDIX A: EVALUATION OF THE SLOPE IN 
THE WEAK MIXING LIMIT 

Here we will describe how to evaluate the slope of the 
GTL (Eq. ([T7}>) at x. = 0. The procedure for x b = 
is the same. The corresponding formulae are obtained 
by interchanging the particle indices b <-> s. Let us fur- 
ther remark that the explicit specialization on a certain 
model system occurs only on the level of the static in- 
put for MCT. Thus, the MCT formulae presented below 
can be directly translated and applied to arbitrary binary 
mixtures such as soft sphere mixtures or binary discs in- 
cluding dipolar interactions 14 1 . 



1. Rewriting the mode coupling functional 

For the following, it is convenient to rewrite the mode 
coupling functional as 



T = nx^ x Tx^ 1 



(Al) 



where the elements of the matrix x are defined by xff = 
x a S a p. As can be read off from Eqs. (J4j)-(l6j> , T has a 
binlincar functional dependence on the matrix c of di- 
rect correlation functions, and shows no further explicit 
dependence on the control parameters. J- can be consid- 
ered as a special case of a more general functional 



F[X,Y]=F[c,c;X,Y], 



(A2) 



^[a,b;X,Y] = £ dp dg 

( 47r r„,^ fl »-/o J\k- P \ 



xv; 



a' ,0' ,ct" ,/3' 
af3;a' f3' ,a" j3" i 



k;p,q 



a,b] 



xX^'Yf^", (A3) 



of a pair of derivatives of the separation parameter of the 
form 9cr/9(A£j)| A |- = jj. Since a follows from a (Eq. p^]) ) 

by linearization around £ c , we can write 



da/d(A^)\ A ^ = da/d^ =ic 



(A8) 



Only, those quantities on the r.h.s of Eq. p6[) without 
the superscript c arc differentiated. Then, all quantities 
in the resulting formula have to be taken at the critical 
point £ c . For the following, we drop the superscript c for 
convenience. With Eqs. (fll)|) , (|A6p and (|A8| we obtain 
explicitly 

da/d(A^)\ Ai=5 = n-^dn/dmHliS-F} 

xnx- 1 T[F,F]x- 1 {S - F}) 
+ (H\{S-F} 

xnx^QlF, F]x~ 1 {S - F}) 
+ (HUS-FjS-^dS/d^} 

xnx^TiF, F}x- 1 {S - F}) 
+ (H\{S-F} 

xnx- 1 T[F,F}x- 1 {dS/d^ i }) 
+ {H\{S-F}n{dx- l /d^} 

xT[F,F]x- 1 {S - F}) 
+ (H\{S - F}nx- l P[F,F] 

xidx^/dtijiS-F}). (A9) 

Note that for a one-component model we have x = x^ 1 = 
1 and thus dx -1 /d^t — 0. Let us further remark that the 
first scalar product on the r.h.s. of Eq. (|A9[) is nothing 
but the well-known exponent parameter A = (H\{S — 
F}F[F,F]{S - F}). 



V, 



a (3; a' j3 \fx" (3" 



k;p,q 



[^b] = ^v^''[a]v^''[b], (A4) 



_ Q/ 3 7f (fc 2 + p 2 - g 2 )zf8 ai + (fc 2 - p 2 + g 2 )z^6 af3 

V kpg l Z \ { Ak 2 p 2 _ ( fc 2 +p 2_ g 2)2](3-d)/4 

(A5) 

Hence, for fixed X and Y and some arbitrary external 
control parameter £j we can write 



(dr/dZi)[x,Y] = g[x,Y], 



(A6) 



G[X, Y] = F[dc/d£i, c; X, Y] + T[c, dc/dfr, X, Y]. 

(A7) 



2. Derivatives of the separation parameter 

a. General case 

For a general model system, the calculation of the slope 
of an arbitrary GTL (Eq. (fT5)0 requires the calculation 



b. Weak mixing limit 

We specialize Eq. (|A9I) to evaluate Eq. p?|) at x s = 0. 
Let us start with summarizing some important properties 
of S, F, T and Q. By definition, for x s — > the elements 
S and F satisfy 



Sf = <D{x 8 ) if (a,/3)^(M), 



(A10) 



Ff = 0{x s ) if (a,(3)^(b,b) 
Due to the Kronecker deltas is Eq. (|A5|) . we also have 



?f[F,F]=0{x B ) if (a,0)^(b,b), 
gf[F,F]=0(x s ) if (a,f3)^(b,b). 



(All) 



For the following, we assume Taylor expansions for n, c, 
S, F, H and H in powers of x 8 around x s — of the 
form 



Z = Z^+x s Z^+0(x 2 s ). 



(A12) 
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Equation (IA10|) implies 

(S(°))f = if (a,/3)^(b,b), 
(FW)f = if (a,f3)^(b,b). 



(A13) 



The Taylor expansions of T[F, F] and Q[F,F] needed 
below read explicitly 

P[F,F] = F (a) [F^\F^]+x s g (0) [F^\F^] 

+ 2x s F i0) [F(°\F^} + O(x 2 s ), (A14) 

Q[F, F] = d (0) [-F (0) , F (0) ] + 0{x s ), (A15) 
where the leading order functionals are given by 



T ia) [X,Y] =F[c(°\c^;X,Y}, 



(A16) 



g (0) [X,Y)= T[ C W , c (°) ;X,Y]+ F[c<V , c« ; X, Y]. 

(A17) 

Equation (|A13p and the Kronecker deltas is Eq. (|A5|) 
imply 

(^ (0) [F(°),F(°)])f ^0 if (a,0)^(b,b), (Alg) 
(e (0) [F(°),F(°)])f = if (a,p)?Q>,b). 

A further important implication is the fact that 
(^[ F (o) )F (i)])«/3 is not dependent on if 
(a,/3)^(b,b). _ 

Now we consider the numerator in Eq. (|17p . It follows 
from Eq. (|A9|) by choosing £j = a; s . Let us focus on the 
scalar product in the first term on the r.h.s. of Eq. (|A9|) . 
The factors {S - F}x-\ n, F[F,F] and x~ 1 {S - F} 
have all well defined limits for x s — > which can be cal- 
culated independently. Hence, the limit of the second 
argument of the considered scalar product also exists. 
Thus, the x s — > limit of the first argument of the scalar 
product, namely that of H , can be performed indepen- 
dently with H^> as result. Because of Eqs. (|A14I) and 
(IA18p . the final result for the x s — > limit of the first 
term on the r.h.s. of Eq. (|A9[) depends only on the ma- 
trix elements with indices (a,/3) = (b,b). We can write 
the result explicitly as {n^ /n^}\^ where 

A (0) = „(0)^(H(°))^{(S(°))^ -{F^) bb } (A19) 

k 

x(P i0) [F(°\FW})f{(SW)f - (F^)f} 

is nothing but the well known exponent parameter of the 
corresponding monodisperse MCT model [5|, [30|. The 
second term on the r.h.s. of Eq. (|A9[) can be discussed 
similarly, here Eqs. (|A15[) and (|A18|) lead to 



Y,(H^)l b {(S^)t b -(F^)t b } (A20) 

k 

{ g {0 \ F V\F(%f{(sW)f-(FM)f}. 



The treatment of the remaining terms in Eq. (|A9[) is 
somewhat more tedious. For this purpose we write 
the matrix products occurring as second arguments of 
the scalar products explicitly in components. By using 
Eqs. (|A10p and (|A1 1|) we realize that all the inverse pow- 
ers of x s stemming from x^ 1 and its derivative with re- 
spect to x s can be compensated by other factors which 
are of the order x s . Hence, the x s — > limits for all 
matrix products occurring as second arguments of the 
scalar products exist. Thus, for each scalar product, the 
x s — > limit of H can be performed independently yield- 
ing H(°\ The final result for numerator in Eq. (jTTJ) eval- 
uated at x s — can be written as 

da/dx s \ Xs=0 = {2 + ? i( 1 Vn (0) }A(°) + M W 

+ (£<®\{A<® (A21) 



(^ (0) )f 



i^>{(snr~(F (0) )t} 



1 (0)/r< ? (0)N 
^-(°)r 



(A22) 



x(^"[F^,F^]) bb (S^)f 
2n(°>{(S<°>)f - (FM)f } 

2n(°){(s( i ))r - (F (1) )n 



x(^ (0) [ir(°),FW])f{(S( )) 
4n(°){(SW)r-(FW)n 
x(t {0) [FM,FW]y k °{(SM)f 
2n^{(S^Tk S -(F( 1) )r} 



{F (0)^ } 

(F (l))f } 



x(F {0) [FW,FM})i s (SW) 



S0 

fc > 



(A23) 



(K^)f = l-(F(°))f/(S(°)) h fe 6 , 

( ^(o ))fcfeS = (f(0))»(5W)k-/(5(°))»-(J'W)^, 



= 0, 



fe ' 



(A24) 



x{(S(°))f-(F(°))f} 

+ 2n(°>(S( 1 ))n^ ( V ( °\* ,(1) ])f 
x{(S(°))f-(F(°))f} 

+ 2n(°>(s«)r(^ ( V (0 \* ,(1) ])r 

x{(S«)f -(F«)f}, (A25) 



Due to the statement below Eq. (|A18|) . the final re- 
sult, Eq. (|MT|) . does not depend on (-F (1) )^ b . The 
term 2A'- *' results from the 6&-elements to the last two 
scalar products in Eq. (|A9|) . The matrix represents 
the contribution of the third term in Eq. (|A9|) where 
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is nothing but the x s -> limit of {S - F}S~ 1 
while is the corresponding limit for the expres- 

sion {dS / dx s \nx~ 1 T[F , F]x~ 1 {S — F}. All remaining 
quantities are summarized to the matrix A^°\ 

Let us now consider the denominator in in Eq. (117|) 
which follows from Eq. (|A9|) by choosing £j = p. Since 
Eqs. (|A10p and (|A11|I remain valid if one replaces the 
corresponding quantities by their derivatives with respect 
to p and since dx~ 1 /dtp = 0, the final result depends 
only on the &&-matrix elements. Thus, the denominator 
in Eq. (fT7|) taken at x s = follows directly from the 
separation parameter of the monodisperse system. It is 
a positive constant. 



3. Slope of a critical line 

The explicit results above allow us to define a proce- 
dure for the calculation of the slope of a GTL at x s = 0. 
It consists of five steps. 



a. Calculation of the critical point 

The first step is the determination of the critical pack- 
ing fraction ip^ and the corresponding NEPs (_F C '(°))^ & 
by using the corresponding one-component model of big 
particles. In the following, all quantities have to be 
taken at ip — ipQ, the critical packing fraction of the 
one-component system. The denominator in Eq. (|17j) 
taken at x s = also follows directly from the separation 
parameter of the monodisperse system. It is a positive 
constant which we calculate by numerical differentiation, 
for simplicity. 



x s c^ + . . . , h analogous and x = x^ + x s x^ into 
Eqs. ([!|]) and (flU|) leads to the equations for c^ n ' and h^ n > 
which have to be solved recursively. For n = and n = 1 , 
they read 



h (o) =c (o) + n (o) c (o) x (o) /l (o) 

with the zeroth order PY closure 

( h (0))a/3( r ) = _ h r <{R a + R ), 
(c(°>)°"(r) = 0, r > (R a + R p ), 

and 



(A28) 



(A29) 



+c (o) x (D h W +c W x ^) h (D } (A3Q) 
with the first order PY closure 



(fc(i))«/J( r ) = o, r < (R a + Rp), 
( c (D)«/3( r ) = 0, r> (R a +Rp). 



(A31) 



Furthermore, we have (x (0) )| b = 1, {x {1) )f = -1, 
^ x (i)^ss _ j anc j & y[ other components are zero, and 
and are given by 



n<°> = (<pd)/(fi«,Jtf), 

„(l) = n(°)(l-(i? s /i? b ) d ). 



(A32) 



Note that Eqs. (|A^9|) . (|MT|l and (|Al2|) are the only ex- 
plicitly model dependent equations. Hence, the proce- 
dure can be easily extended for both to arbitrary binary 
mixtures and to closure relations different from PY. Let 
us further remark that (c^) h £ and are nothing 

but the direct and total correlations functions for the 
one-component system of big particles. 



b. Calculation of the static structure 

and entering into da/dx s \ Xe= o trough 

Eqs. (|A19[) - (|A25[) can be easily determined from and 
by using Eq. (J7J). The result reads 

(S< ))? = l/{l-n< >(c<°>)?}, 

(S^)l s = 0, (A26) 

(s (0) )r = 0, 



(s«)r = i, 

{S (l) )bk S = n ( ) {S W)f(cW)f, 

(SW)f = {(S^)l b } 2 

x{nW[(cW)f + (c<°>)jfl 

_( n (0)_ n (l))( c (O))66 

-(« (0) ) 2 [(c( o ))^( C (°))r-(c (0) )r(c (0) )f]} 
- (s(°))! b {i+™(°)( C (°))r}. 

(A27) 

Hence, in the second step we have to determine c' ' and 
cW. Substituting n = + x s n^ + c = + 



c. Calculation of the critical nonergodicity parameters 

Beside {F^)f , the evaluation Eqs. ([A"l9l ) -(TAl5 ]) re- 
quires also (F^)^f and (F^)l s as input. It is straight- 
forward to derive the equations for these quantities from 
the fixed point equation F = T[F] following from 
Eq. (jTTj) by considering the limit x s — > 0. We obtain 

= l-{l + 2n(°)(^ (0) [F(° ) ,F( 1) ])^}- 1 , (A33) 
(F«) h fe s = 2n(°>(S«)^ ( V ( °\* ,(1) ])r 

x{i-(F«)n 

+ 2n<°>(S< >)^(:F ( V ( °\* ,(1) ])fc S 
x{l-(F«)n 

+ n(°)(5( ))f(^ (0) [F(°),F( )])^ 

x{(S( 1 ))^-(FW)t 8 }- (A34) 

Since {F^) b £ have already been determined in the first 
step, Eq. (|A"33)l allows to calculate {F^)% s . The r.h.s. 
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of Eq. (|A33[) does neither depend on {F^')^f nor on 
(F (1) )bb. The are nothing but the tagged par- 

ticle NEPs for a single small particle in the fluid of the 
big particles. Finally, Eq. (|A34[) allows us to calculate 
(F^ 1 ')^, since it is not dependent on (F^)^ due to the 
statement below Eq. (|A18|I . 



Now, and the corresponding adjoint map C7(°) allow 
us to calculate the eigenvectors and obeying 

the normalization, 



k 



(A38) 



d. Calculation of the critical eigenvectors 



The evaluation Eqs. (|A19[) - <|A25ji requires the zeroth 
order left eigenvector as last input. For its unique 

determination, also the zeroth order right eigenvector 
J?( 0) is needed. For x s -> 0, Eq. (TjTJ]) reduces to 

C<®[Y] = 2n^M^T i0) [F ( -°\Y}N ( - \ (A35) 



E 

fc 



{(HW)ff/{(SW)f - (F<°>)?} = 1. 



(A39) 

While for only the 66-elements are nonvanishing, 

if has nontrivial contributions for all particle indices. 
(J3-(o))W> an d (H^)f are the eigenvectors for the one- 
component model of big particles. 



(M(°>)S* 
(M(o))fe 

(M(°))f 
(M(o) )r 



( 5 (0)' 

0, 

i-(f«; 



\6b 



(J .(0) ) » j 



(A36) 



it ' 



Calculation of the slope 



(jV(°))j* 
(7V(o))6 S 

( Ar(o) )r 



Nfcb 



(_p(o))fcfc 



0, 

1 - (FW 



fe ' 



\s6 



(A37) 



Now, we have determined all quantities for the evalu- 
ation of Eqs. (|A19[) - (|A25|) and are able to calculate the 
slope of the GTL by using Eqs. QH) and l[A"8j) , 



)k ■ 
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